Spectroscopy on two coupled flux qubits 
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We have performed spectroscopy measurements on two coupled flux qubits. The qubits are 
coupled inductively, which results in a erf erf interaction. By applying microwave radiation, we 
observe resonances due to transitions from the ground state to the first two excited states. From 
the position of these resonances as a function of the magnetic field applied we observe the coupling 
of the qubits. The coupling strength agrees well with calculations of the mutual inductance. 

PACS numbers: 03.67.Lx, 74.50.+r, 85.25.Cp 



I 

-a 
c 

o 
o 



> 

(N 
00 

o 
m 
o 



I 

C 

O 

o 



X 



Quantum computers manipulate quantum information 
contained in the states of interacting two-level systems 
called quantum bits or qubits. Quantum gates require 

!|ubit-qubit coupling that preserves quantum coherence 
lj. Qubits have been implemented in various systems. 
However, the requirement for upscaling makes solid state 
implementations highly attractive. Several single solid 
state qubits have been realized using superconducting 
Josephson junction circuits 0, IE 0, la, IE ■ As a next 
step coupled multiple qubits need to be studied. Re- 
cently charge quantum dynamics in two coupled Joseph- 
son qubits has been observed Q . Spectroscopic measure- 
ments on coupled charge qubits have been performed as 
well @. 

In this letter we present spectroscopy measurements 
on two coupled flux qubits |2j. Our flux qubit consists 
of a small superconducting loop interrupted by three 
Josephson junctions with a small junction capacitance C . 
When approximately half a flux quantum $o is applied to 
the ring the two classical states correspond to clockwise 
and anti-clockwise circulating currents. The circulating 
current of the qubit generates a magnetic field into |J.) 
or out of |t) plane. Due to the large charging energy 
Ec = e 2 /2C the barrier between these classical states 
is low and tunneling between the states will take place. 
When exactly half a flux quantum is applied, the two 
eigenstates are symmetric and antisymmetric superposi- 
tions of the two classical states. The system is described 
by the effective Hamiltonian H = ha z + ta x where a z,x 
are Pauli spin matrices and h — I p (& — l/2<f>o), with 
I p the circulating current and $ the flux applied to the 
qubit. The tunneling amplitude t depends exponentially 
on the ratio Ej/E c Hfl. Here Ej = I <& /2tt is the 
Josephson energy, with I the critical current of one junc- 
tion. Because the basic states of this qubit are flux states 
it is very insensitive to charge noise. The energy level 
repulsion due to the tunneling amplitude t was observed 
with spectroscopy 0. Recently also coherent oscillations 
of such a flux qubit have been observed . 

The easiest way to realize inter-qubit coupling for flux 
qubits is by magnetic interaction. This coupling results 




FIG. 1: Scanning microscope image of the two coupled qubits 
surrounded by the DC-SQUID. A part of the microwave cou- 
pler is visible on the right. The width of the qubits is 5 fj,m. 



from the magnetic flux that the circulating current in 
one qubit induces in the second qubit. The interaction is 
of Ising form <j\ erf, since the clockwise |J.) and the anti- 
clockwise \\) circulating states give a z-component field. 
The coupling strength j is given by 



j = 2MI pA I Pi2 



(1) 



where M is the mutual inductance between the qubits. 
MI Pt i is the generated flux in qubit 2 by the current 
in qubit 1 leading to an increase in energy of MI p ^I p p 
in the second qubit. The sign of the coupling is such 
that it favors anti-parallel qubit states. In the two level 
description the self induced flux can be ignored because it 
adds a constant energy term. Figure Q shows a scanning 
electron microscope image of the two qubits. In order to 
enhance the qubit-qubit mutual inductance M the two 
qubits share one branch. 

The two qubits are now described by the total Hamil- 
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FIG. 2: Energy levels of two identical coupled qubits as a 
function of the externally applied flux. Far away from the 
degeneracy point \h\ 3> t, the energy states resemble classical 
states. The coupling j shifts the energy of the ferromagnetic 
states ITT), up, while the energy of the antiferromagnetic 
states |t|), ||t) is lowered. In the vicinity of the degener- 
acy point \h\ < t, the classical states mix due to the tunnel 
coupling t. Due to the symmetry of the problem, the anti- 
symmetric energy state |||) — ||J.) does not mix with sym- 
metric states. 
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Here hi ~ I p ,i($i ~ 1 /2<3? ) and ti is the tunneling 
amplitude of the first qubit. Similar expressions hold for 
the second qubit. 

The energy levels and energy eigenstates of two cou- 
pled qubits can be calculated by diagonalizing the Hamil- 
tonian J2J. A discussion for coupled solid-state qubits is 
given in [13J . Figure [21 shows the energy levels for two 
identical qubits. The parameters are I p .\ = I p .2 = 1 
GHz / m$o a nd j = h = ti= 1 GHz. Far away from the 
degeneracy point (\h\ ^> t), the energy states closely re- 
semble the classical states |J.T), ||4), | T T ) - The two 
antiferromagnetic states |||), ||t) (with opposite circu- 
lating currents in the qubits) are degenerate. The cou- 
pling j reduces the energy of the antiferromagnetic states 
and increases the energy of the ferromagnetic states ||t)) 
|4I). In the vicinity of the degeneracy point (\h\ < t), the 
states are superpositions of the classical states due to the 
tunnel coupling t. Due to the symmetry of the problem 
the anti-symmetric energy state |||) — ||J,) does not mix 
with the other states and its energy is independent of the 
magnetic field applied. 

In practice the two qubits have different parameter val- 
ues due to fabrication limitations. Figure EKa) shows 
the energy levels for two unequal qubits. The param- 
eter values are those of the actual sample (as obtained 
below) . Because the qubits have different circulating cur- 



FIG. 3: (a) Energy levels of two different coupled qubits as a 
function of the externally applied flux. Due to the difference 
between the two qubits, the degeneracy of the two antiferro- 
magnetic states HI), HI) is lifted 

(b) Transition matrix elements |Ti/| = \ f\ af + a 2 |^i) | 
for the transitions from the ground state to the three excited 
states. The matrix element for the transition from the ground 
state to the third excited state is very small except for the 
vicinity of the degeneracy point. The parameter values are 
taken from the values obtained by fitting the experimental 
data (Fig. 0J. 



rents / Pi i ^ / p ,2 the degeneracy of the antiferromagnetic 
states is lifted. 

The signal of the two qubits is measured with a DC 
superconducting quantum interference device (SQUID). 
The SQUID surrounds the two qubits (Fig. ^) and mea- 
sures the magnetization, i.e. generated flux, of the two 
qubits together. The magnetization of the two qubits is 
given by the slope of the energy levels as a function of 
the magnetic flux (Fig. [3 EJ. On the right-hand side 
of figure ^a) part of the microwave coupler is visible. 
When the frequency of the microwave matches the en- 
ergy difference between two energy states of the coupled 
qubits the radiation induces transitions between the two 
states. The microwave coupler induces approximately 
the same amount of flux in the two qubits. Therefore the 
transition matrix element is of the form a \ + erf . Figure 
|3tb) shows the absolute value of the transition elements 
T if = (9f\af + 

Measurements were performed in a dilution refriger- 
ator at a base temperature of 20 mK. We measured 
the switching current of the DC-SQUID. This was done 
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by ramping the bias current and recording the current, 
where the SQUID switched to a voltage state We re- 
peated this procedure and averaged the switching current 
over typically 1000 measurements. We varied the field in 
a 20 m$ range centered at 1/2 $ . The SQUID signal 
shows a characteristic step (Fig. Ufa)). The ground state 
of the two qubits changes from a state via a super- 
position of all states at 1/2 $ to the |J,|) state. Applying 
microwaves leads to dips on the left of the step and peaks 
on the right of the step (Fig. EJa)). At these positions the 
energy of the microwaves matches the difference between 
energy levels and induces transitions between them. As 
the system of the two qubits is excited continuously an 
incoherent mixture is formed and the SQUID signal is 
the average of the magnetization of the two levels. We 
have studied the positions of these resonances as a func- 
tion of the applied microwave frequency. Figure EJb) 
shows the frequency of the resonance versus the position. 
Far away from 1/2 $ the peaks and dips approximately 
follow straight lines. Here the resonances correspond to 
transitions between states where one of the two qubits 
is flipped. One observes that the slopes of the two res- 
onances and consequently the persistent currents of the 
two qubits are different. The dotted lines inEJb are lin- 
ear fits to the data above 8 GHz for one of the qubits. 
From these lines an estimate can be made of the coupling 
strength j. As the ground state is ferromagnetic and has 
its levels shifted up by j, while the excited state is anti- 
ferromagnetic with levels shifted down by j, the intersec- 
tion should occur at a level of —2j. For the dotted lines 
one finds a crossing at a level of -0.98 GHz, indicating a 
value j — 0.49 GHz. The transitions for the other qubit 
yield a less accurate estimate as can be understood later 
from the full parameter fit (the higher value of ti leads 
to strong rounding). One observes that the resonance 
peaks and dips cross at two points at opposite sides of 
the degeneracy point 1/2 $o with the crossing points at 
~ -1.5 m$o for the peaks and ~ +2.2 m$o f° r the dips. 
This is a peculiar feature of the quantum nature of the 
two qubits. If the qubits were classical (t\ = t% = 0) the 
energy levels in figure IHIa) would follow straight lines. 
The first and second excited state lines thus could cross 
at one point only. 

In order to obtain all coupled qubit parameters we 
did a least-squares fit to the difference of energy lev- 
els calculated with the coupled qubit Hamiltonian J2J. 
Here we make the assumption that the magnetic field 
is homogeneous for the qubits and the SQUID. The 
least square fit involves six fitting parameters: the per- 
sistent currents and I p ^, the tunneling amplitudes 
t\ and ti-, the inter-qubit coupling strength j, and the 
relative difference in area of the two qubits a. From 
the fit we obtain the following parameters for the two 
qubits: I pA = 512nA ± 6nA, t x = 0.45 GHz ± 0.2 GHz 
and I Pt2 = 392 nA ± 5nA, t 2 = 1.9 GHz ± 0.1 GHz. 
The persistent current is directly dependent on the crit- 
ical current of the Josephson junctions 10]. The dif- 
ference in persistent currents of the two qubits is due 
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FIG. 4: Spectroscopy measurements, (a) SQUID signal, i.e. 
magnetization of the two qubits, as a function of the applied 
magnetic field. The microwaves induce transition from the 
ground state to the first and second excited state, which ap- 
pear as dips and peaks in the magnetization. The outer peaks 
and dips are due to the first excited state and the inner due 
to the second, (b) Microwave frequency versus peak and dip 
position. In the vicinity of the degeneracy point ($ ~ l/2$o), 
the resonances deviate from straight lines because the energy 
levels are formed by superposition of the classical states. The 
straight line is a fit of the two qubit Hamiltonian (Eq. |5J to 
the data. The inset compares the data obtained in the vicin- 
ity of 3/2$o with the data obtained at l/2$o- The shift is 
due to the influence of the different qubit areas which is three 
times higher. 



to the spread in critical currents in the fabrication pro- 
cess. While the tunneling amplitude t depends expo- 
nentially on Ej/Eq, the typical spread in the junction 
parameters can easily lead to a difference as large as ob- 
served. The qubit coupling strength obtained from the 
fit is j = 0.50 GHz ± 0.03 GHz and agrees well with re- 
sult from the analysis above. The difference in qubit area 
is cr = 0.027% ± 0.004%. If the areas would differ more 
than a few parts in a thousand, the degeneracy points 
of the two qubits would be so far apart that one would 
measure a double step. 

The coupling strength j can also be evaluated using 
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Eq. Q). From the geometry the mutual inductance can 
be calculated, yielding M =0.86 pH. Using the values 
persistent currents obtained in the fitting procedure, we 
get a coupling strength of j = 0.52 GHz, which agrees 
well with the result above. 

We have also performed the spectroscopy measure- 
ments in the vicinity of 3/2 $o- At this point everything 
is equivalent to the situation at 1/2 $o except for the 
flux bias due to the difference of the areas. The effect of 
the different qubit areas on the energy levels is now three 
times higher. The inset in figure^b) compares the data 
at 3/2 and 1/2 <£>o at high frequencies. The peaks due to 
the first excited state are shifted to lower magnetic field 
values and the peaks of the second excited state to higher 
values. From the magnitude of the shift one can calculate 
the difference of the qubit areas to be a = 0.03%, which 
agrees with the result from the fit above. 

We have only observed transitions from the ground 
state to the first two excited states, but not to the third 
excited state. Figure GJb) shows that the transition am- 
plitude T 3 for this resonance is very small except for the 
vicinity of the degeneracy point 1/2 $o. However, near 
the degeneracy point the ground state and the third ex- 
cited state have similar slope and therefore similar mag- 
netization. Consequently we are not able to observe the 
transition with the SQUID. Figure |3fb) also shows that 



the transition amplitude to the first excited state is larger 
than to the second excited state. In agreement with that 
prediction, one observes that the outer resonance peaks 
and dips are broader than the inner ones (Fig. Ufa)). 

The previous analysis makes clear that the spectro- 
scopic data are fully consistent with the two-qubit Hamil- 
tonian of Eq.J2J. This Hamiltonian, in turn, opens the 
possibility of well-chosen one and two-qubit operations 
that lead to controlled entanglement. The new results 
support the notion that superconducting flux qubits can 
be used to study entanglement in macroscopic quantum 
systems and for the development of non-trivial two-qubit 
gates such as the controlled- not. 

In summary, we have performed spectroscopy mea- 
surements on two coupled flux qubits. The mutual in- 
ductance between the two qubits leads to Ising coupling 
uf erf. The observed resonances agree very well with the 
two qubit Hamiltonian (Eq. The magnitude of the 

qubit coupling agrees well with the calculation of the mu- 
tual inductance between the qubits. 
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